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1. Introduction 

The predictive power of QCD for high-energy hadronic scattering observables rests upon 
our ability to compute the corresponding partonic cross sections in perturbation theory. 
Factorization theorems [1] assert that for many such cross sections mass divergences may 
be subtracted, or mass-factorized, in a process-independent way, with any additional finite 
subtraction constants defining the mass-factorization scheme. 

In many cases, the predictive power of such perturbative series is imperiled by the 
systematic occurrence of finite but large terms at higher orders. Resummation attempts 
to restore predictive power by organizing classes of large terms to all orders, leaving a per- 
turbative series for the remainder with much better convergence properties [2]. For cross 
sections, the large terms almost always take the form of logarithms of ratios of kinematical 
scales. In particular, threshold resummations [3, 4] organize logarithmic enhancements 
singular at partonic threshold, resulting from imperfect cancellations between real and 
purely virtual contributions to the cross section. As partonic threshold is approached, 
these enhancements are parametrically guaranteed to increasingly dominate the perturba- 
tive contributions to the cross section. Often, however, also constant terms, which do not 
depend on scale ratios vanishing at threshold and which arise predominantly from purely 
virtual diagrams, are numerically important in the cross section. Some of these large con- 
stants originate from the same infrared singularities that give rise to the large logarithms, 
and consequently are resummable. 

Threshold resummations can be derived for partonic cross sections using the procedure 
of refactorization: the Mellin transform of the cross section is expressed near threshold as a 
product of well-defined functions, each organizing a class of infrared and collinear enhance- 
ments. The refactorizations are valid up to corrections which are nonsingular at threshold, 
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and thus suppressed by powers of the Mellin variable N at large N. Terms independent of 
TV can then be treated, in principle, by the same methods used to resum terms enhanced by 
logarithms of N. For the Drell-Yan cross section and the deep-inelastic structure function 
F2 such a refactorization was achieved in Ref. [3] . The resulting resummation of constant 
terms established to all orders the earlier observation [5] that, in the DIS-scheme Drell-Yan 
cross section, the largest constants are related to the ratio of the timelike Sudakov form 
factor to its spacelike continuation. By solving an appropriate evolution equation [6], an 
exponential representation for this form factor, and thus also for the ratio, was derived [7], 
generalizing the observation of [5] to a full nonabelian exponentiation. Comparison with 
exact two-loop results [8] showed in that case that AMndependent contributions at two 
loops are indeed dominated by the exponentiation of the one-loop result, combined with 
running coupling effects [9, 10]. 

In this paper we refine this analysis for the DIS scheme, and we extend it to the 
MS scheme. The challenge in the latter case lies in the fact that the finite subtraction 
constants of this scheme are not related to a physical scattering process involving the 
electromagnetic quark coupling at lowest order. Therefore, ratios of form factors do not 
naturally occur in the MS scheme. The practical prevalence and relative simplicity of 
the MS scheme would, however, make such an organization desirable. We will show that 
the refactorization formalism of Ref. [3] leads to the exponentiation of all AMndependent 
contributions to the inclusive Drell-Yan cross section, both in the MS scheme and in the 
DIS scheme. As a corollary, one may note that all constant terms in the MS -scheme non- 
singlet deep-inelastic structure function F2 have been organized into an exponential form 
as well. Furthermore, it is possible to organize the factorization procedure so that real and 
virtual contributions are individually made finite; one can then disentangle various sources 
of constants, such as 7r 2 terms arising from unitarity cuts and similar terms arising from 
expansions of phase-space related T functions. 

One might object that there is no kinematic limit in which AMndependent terms 
dominate parametrically, so that an organization of such terms cannot be of much practical 
use. In general, our view is that whenever all-order information is available one should 
make use of it, at least to gauge the potential impact of generic higher order corrections on 
the cross section at hand. One should bear in mind that the pattern of exponentiation, even 
for AMndependent terms, is highly nontrivial, and includes all-order information arising 
from renormalization group evolution and the requirements of factorization; for example, 
a considerable fraction of nonabelian effects arising in the Sudakov form factor at two 
loops can be shown to follow from running coupling effects implemented on the (essentially 
abelian) one-loop result. In the present case one could be bolder and argue that, since 
all constants have been shown to exponentiate, using the exponentiated expression should 
provide a better approximation to the exact answer. This is in fact the case for generic 
values of higher order perturbative coefficents, even for asymptotic series such as those 
arising in QCD. It cannot, however, be proven for any particular cross section, although it 
works in practice for the cases that have been tested. At the very least, differences between 
results for the physical cross section with or without the exponentiation of constant terms 
can provide nontrivial estimates of errors due to (uncalculated) higher order corrections. 
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Notice in passing that constant terms are not affected by the Landau pole and thus factor 
out of the inverse Mellin transform needed to construct the physical cross section. 

This paper is organized as follows. In the next Section we review the resummation 
of the quark Sudakov form factor and its embedding in the DIS-scheme Drell-Yan cross 
section. We extend the exponentiation in this scheme to include all constant terms. In 
Section 3 we derive our result for the MS cross section. In Section 4 we compare the results 
of the exponentiations to finite order calculations, while Section 5 contains our conclusions. 
In an Appendix we present renormalization group studies of certain auxiliary functions. 

2. Exponentiation in the DIS scheme 

Consider the N—th moment of the partonic Drell-Yan cross section, taken with respect to 
z = Q 2 /s, with Q the measured invariant mass and s the partonic invariant mass squared. 
Mass factorization of this quantity, in the DIS scheme, is performed by simply dividing its 
dimensionally regularized (d = 4 — 2e), unsubtracted version by the square of the iV-th 
moment (taken with respect to the partonic Bjorken-x variable) of the non-singlet partonic 
deep-inelastic structure function F2, 

CdisW = . (2.1) 

[F 2 (N,e)} 2 

While numerator and denominator are each infrared and collinear divergent, their ratio 
is finite to all perturbative orders [1], so the e dependence of the left hand side can be 
neglected. Dependence on the hard scale Q and on the strong coupling a s (/U 2 ) will generally 
be understood. 

Let us begin by identifying how the Sudakov form factor arises in uj(N, e), before mass 
factorization, following the reasoning of Ref. [3]. One observes that near threshold (i.e. at 
large N) u(N, e) can be refactorized according to 

lu(N, e) = \H DY \ 2 if>(N, ef U(N) + 0(1/N) . (2.2) 

The ip{N,e) and U(N) functions contain the ./V dependence associated with initial state 
radiation at fixed energy and coherent soft radiation, respectively. They are well-defined 
as operator matrix elements, as given in Ref. [3], and are calculable in perturbation theory. 
Both functions are gauge-dependent, but their product in Eq. (2.2) is not; implicitly, we 
will be working in axial gauge. Using gauge invariance and renormalization group (RG) 
arguments, one can show that both the parton distribution ip(N, e) and the eikonal function 
U (N) obey evolution equations which can be solved near threshold in an exponential form, 
up to corrections suppressed by powers of TV. The function \Hdy\ 2 , collecting all hard-gluon 
corrections, has no iV dependence and may be determined by matching to exact calculations 
order by order. Divergences are only present in the parton distribution function ifj(N,e). 

To identify the Sudakov form factor, it is useful [3] to separate virtual (V) and real 
(R) contributions to the resummed ijj and U functions, according to 

Y>(7V,e) = K(e)il> R (N,e) 
U(N) = U v (e)U R (N,e) , (2.3) 
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where 11(e) is the real part of the residue of the quark two-point function in an axial gauge. 
Note that this further factorization makes sense only because the functions ip and U are 
known to exponentiate to the present accuracy. Using the analysis of Section 8 of Ref. [3] 
we can now write 

u(N,e) = \H DY 1Z(e)^/U v (e)\ 2 ^ R (N, e) 2 U R (N) + 0(1/N) 

= \T(Q 2 ,e)\ 2 ij R (N,e) 2 U R (N,e) + 0(l/N) , (2.4) 

so that all virtual contributions are expressed in terms of the electromagnetic quark form 
factor T(Q 2 ,e). In fact, the residue of the quark propagator coincides with the virtual jet 
function summarizing virtual collinear contributions to the form factor, while the square 
root of the virtual eikonal function appearing in the cross section is responsible for the soft 
enhancements of T(Q 2 ,e). We have thus indentified the dimensionally regularized time- 
like Sudakov form factor in the refactorized, unsubtracted Drell-Yan cross section. Near 
threshold, the only remaining contributions to the cross section come from real radiation, 
and are summarized by the real parts of the tp and U functions. At this point one can 
already observe that ujdy(N, e) exponentiates up to corrections suppressed by powers of N: 
the exponentiation of the form factor in dimensional regularization was proven in Ref. [7], 
while the exponentiation of ij) R and U R to this accuracy was proven in Ref. [3]. Specifically, 
the real part of the fixed-energy parton density ip R (N, e) can be written as 

4r(N, e) = exp dz f-^- J 1 H (a ((1 - y) 2 Q 2 ) , e) j . (2.5) 

Similarly 

U R (N, e) = exp j - J dz ^— 9u (a ((1 - z) 2 Q 2 ) , e) j . (2.6) 

Note that in writing Eq. (2.5) and Eq. (2.6) one makes use of the fact that both functions 
are renormalization-group invariant, since real emission diagrams do not have in this case 
ultraviolet divergences. For the function ip R this is a consequence of the fact that it fixes the 
energy of the final state, so that transverse momentum is also limited; for the function U R 
it is a consequence of the structure of nonabelian exponentiation, as discussed in Ref. [3]. 
The functions and gjj are both finite at one loop in the limit e — > 0: all IR and collinear 
singularities are generated by the integrations over the scale of the <i-dimensional running 
coupling, which at one loop is given by 

" (^2' as<J " 2 )' e ) =a s(f j2 ) 

and will often be abbreviated by a(£ 2 ), as in Eqs. (2.5) and (2.6). 

A similar refactorization can be performed on the deep inelastic structure function F 2 
[3]. One finds 

F 2 (N, e) = \H D1S \ 2 X (N, e) V(N) J(N) + 0(l/N) . (2.8) 

Here the parton distribution x(iV, e) has the same soft and collinear singularities as the 
one adopted for the Drell-Yan process, but, according to the general strategy of Ref. [3], it 
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fixes a different component of the incoming parton momentum (the plus-component) and 
is computed in a different axial gauge. It exponentiates in a form similar to Eq. (2.5), but 
with a different function k x appearing in the exponent. Note that and k x must (and do) 
differ at one loop only by terms of order e 2 , so that the ratio tpR/XR ma Y remain finite. The 
soft function V(N) summarizes the effects of coherent soft radiation in the DIS process, 
and exponentiates in a form similar to Eq. (2.6). Finally, J(N) contains the effects of final 
state collinear radiation emitted by the struck parton, and separately exponentiates in a 
somewhat more elaborate form which will not be needed here. Separating real and virtual 
contributions as above we find 

F 2 (N,e) = \H ms \ 2 \K(e)V v (e)\ X R(N,e)V R (N,e)J(N)+0(l/N) 

= \H ms \\T(-Q 2 ,e)\ X R(N,e) y/W^)V R (N,e) J(N) + 0(1/N) . (2.9) 

To this extent virtual contributions to F 2 have been organized already in [3]. We now see 
that all virtual contributions can be organized in terms of T(— Q 2 , e), by observing that the 
purely virtual part of the light-cone distribution x is identical to the purely virtual part 
of the outgoing jet J. Both consist essentially of the full two-point function for a lightlike 
fermion. Note that both virtual jets in Eq. (2.9) are computed with the same gauge choice. 
Gathering all virtual parts, one finds then 

F 2 (N,e) = \T(-Q 2 ,e)\ 2 XR (.N,e)V R (N,e) J R (N,e) + 0(1/N) . (2.10) 

Again, exploiting the results of Refs. [3, 7] this form of the refactorization is sufficient to 
prove the exponentiation of the full cross section up to corrections suppressed by powers 
of N. In fact F 2 now involves, to this accuracy, only the form factor, and a product of 
real functions which have been shown to exponentiate by using their respective evolution 
equations. 

This result can be further verified in the following way. A comparison of Eqs. (2.9) 
and (2.10) implies that -ffois itself acquires an exponential form. In fact, an analysis along 
the lines of [6] reveals that it may be expressed as 

#dis(Q 2 ) = Zh(<*» e) ^ r( ~^ ,e) > (2- 11 ) 

5W (e) 4 0) (QV) 

where = S Z$ and = G 2 Z q are the unrenormalized, dimensionally regularized 
virtual soft function and virtual quark jet function appearing in the factorization of the 
form factor, and where 

Z H = Z S Z 2 (2.12) 

is the UV counterterm function for -ffois- The Sudakov form factor V does not require 
a separate Z factor to cancel QCD UV divergences, by virtue of electromagnetic current 
conservation. Now, each factor in Eq. (2.11) has an exponential form. For the virtual soft 
function this was shown by Gatheral, Frenkel and Taylor in [11, 12]. The unrenormalized 
virtual jet function G^{Q 2 ,() obeys an evolution equation of the same form as the one 
used for the full form factor [6] , which can be explicitly solved in dimensional regularization 
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by the same methods, using as initial condition the fact that all radiative corrections vanish 
at Q 2 = 0. G^(Q 2 ,e) must then exponentiate by itself. Finally, any Z factor arising in 
multiplicative renormalization may be represented in a minimal scheme in terms of the 
associated anomalous dimension 7 = (l/2)d(ln Z)/d\a.fjL as 



cxp ■ 



,Q 2 ^2 



2 7<(*(f 2 )) 



(2.13) 



where again UV poles are generated by integration over the scale of the ci-dimensional 
coupling. 

Turning to the evaluation of Eq. (2.1), we observe that it requires the ratio of Eq. (2.4) 
and the square of either Eq. (2.9) or (2.10). In practice, the expression (2.9) is more 
convenient than Eq. (2.10), because the resulting form (2.15) for u5dis(-W) is a product of 
finite functions. Had one used instead the result in Eq. (2.10), the resulting expression for 
wdis(-W) would have involved cancelling divergences between the real and virtual parts. 
Using then Eq. (2.9), and the additional information that Uy(e) = Vy(e), both being given 
by pure counterterms to the same eikonal vertex, we can write 



wdis(-W) 



1 



\h dis \< 



r(Q 2 ,e) 



r(-Q 2 ,e) 



\Xr(N,€) 



U(N) 



1 



V 2 (N) J 2 (AT) 



(2.14) 



The exponentiation and RG running of the various factors of Eq. (2.14) are described 
in detail in [3], with the exception of the ratio i/jr/xRi which there was exponentiated 
according to Eq. (2.5), but evaluated only at leading order. Running coupling effects on 
this ratio are briefly discussed in an Appendix: they generate a contribution at NNL log 
level at two loops, as well as further ^-independent terms. Furthermore, we are now in a 
position to exponentiate also the one-loop contribution to the matching function, setting 
Hdis = exp (— (XsCf/tt)- Gathering all factors, and formulating the answer according to 
standard notation [4], our result for the hard partonic Drell-Yan cross section in the DIS 
scheme takes the form 




(l-z) 2 < 

(i-z)Q 2 



r 2 A(a s (f)) 



2B (a s ((1 - z)Q 2 )) + D (a s ((1 - z) 2 Q 2 )) 



+ 0(1/N). (2.15) 



Eq. (2.15) resums all terms in the perturbative expansion which contain enhancements of 
the form a n s log fc N, with k < 2n, provided the perturbative expansions of the functions 
A, B and D are known to the desired order in the strong coupling. The perturbative 
coefficients of the functions A, B, D are in fact all known up to two loops [13, 14, 15] 1 . 



1 Concerning the three- loop coefficient A^ 3 \ the n/ -dependent part is known exactly [15, 16], while good 
numerical estimates exist for the n/ -independent term [17]. 
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k=l 



For example, expanding the functions involved as 

oo 

/(*.) = E f (k) 

one needs 

A« = Cf 



= - 
2 



CaCf 



7t) ' 



--C(2) 
18 sv ' 



(2.16) 



n f° F ( 9 



= 



(2.17) 



for resummation to next -to-leading logarithmic accuracy. 

Eq. (2.15) also exponentiates N— independent terms, which have three sources: they 
come from unitarity cuts, as in the analytic continuation of the form factor; or from phase 
space integrations, since for example the parton distributions tp and x differ slightly in their 
phase space measure; finally they arise from the Mellin transformation in the exponent, 
which generates not only logarithms of N, but also contributions proportional to 7e and 
to C( n )- Let us examine in more detail the first two classes of exponentiated constants. 

As far as the ratio of form factors is concerned, one may use the result of Ref. [7] showing 
that the absolute value of the ratio is finite to all orders and exponentiates. To illustrate 
these results, note that the (timelike) Sudakov form factor T(Q 2 , e) for the electromagnetic 
coupling of a massless quark of charge e q is defined via 

r^ Pl ,p 2 ;fi 2 ,e) = (0|J„(0)|pi,p2) = -iee q v(p 2 )^u( Pl ) T (Q 2 , e) , (2.18) 

with Q 2 = (pi + p 2 ) 2 . Based on Refs. [18, 19, 20], it was shown in Ref. [7] that the 
dimensionally regularized Sudakov form factor may be written as an exponential of integrals 
over functions only of a s , 



T(Q 2 , e) = exp 




1 fM 2 wx 2 

K(a s ,e) + G(a(e),e) + ^ (a (X 2 )) 



(2.19) 



The function K in Eq. (2.19) is defined to consist of counterterms only, while G is finite 
for e — > 0. Double logarithms of the hard scale Q arise from the double integral over the 
anomalous dimension 7x(o ; s)- Note that within the framework of dimensional regulariza- 
tion all integrals in Eq. (2.19) can be explicitly performed using the d-dimensional running 
coupling [21], so that the form factor can be expressed in terms of RG-invariant analytic 
functions of a s and e to any order. 

Using Eq. (2.19) one can derive a particularly simple expression for the absolute value 
of the ratio appearing in Eq. (2.15). One finds 

r(Q 2 ,e) 



T(-Q 2 



exp 



G a 



-77- I JBr\1 



cr 



^Q 2 



Performing the scale integrals, at the two-loop level this yields 



T(Q 2 ,e) 



r(-QV) 



(i) 



= 1 + 



a s (Q) 3C(2) 7 ^ 



7T 



(2.20) 



(2.21) 



+ 



a,(Q) 



;C 2 (2)(7^) + |C(2)6 G«(0) + |C(2)71? ) 
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where 



7? = 2C F , 

G {1) (e) = C F Q - | [C(2) - 8] + e 2 [s - |c(2) - |c( 
7 g) = ^C J? (g-C(2))-n / C F (§) =2AW 




(2.22) 



while 60 = (HC^ — 2n/)/3. The anomalous dimension 7^ is the "cusp" anomalous dimen- 
sion of a Wilson line in the MS renormalization scheme [22, 23, 24]. 

Eq. (2.21) illustrates the potential relevance of exponentiation of iV-independent terms: 
first of all, the two-loop contribution is numerically dominated by one-loop effects, both 
through exponentiation and the running of the coupling (the first two terms of the two- 
loop coefficient numerically make up roughly three quarters of the total); furthermore, for 

(2) 

this particular ratio, genuine two-loop effects are given only in terms of "f K , thus they are 
UV-dominated and much simpler to calculate than the full form factor. 

Finally, the function F£>js(a s ) collects constant terms arising from phase space inte- 
grations in the various functions involved in the factorization, as well as from the exponen- 
tiation of the matching function Hp,is- One finds at one loop 



while at two loops some terms can be predicted by taking into account the running of the 
coupling, which yields 



These two-loop contributions should not be taken too literally since, as indicated in 

(2) 

Eq. (2.24), there is at this level an uncalculated contribution SF^j S arising from a pure 
two-loop calculation, which could easily overwhelm the effects which have been included. 
Further discussion of the impact of these two-loop effects is given in Section 4. 

3. Exponentiation in the MS scheme 

As remarked in the introduction, one should not expect that the constants associated with 
the Sudakov form factor in the MS -scheme Drell-Yan cross section can be organized as in 
the previous Section, in terms of a simple ratio of the timelike to spacelike versions of the 
same form factor. The reason, of course, is that the MS quark distribution function is not 
directly related to a physical process involving quark electromagnetic scattering at lowest 
order. We will show that it is nevertheless possible to organize these constants in a closely 
related manner. 

Mass factorization of the Drell-Yan cross section in the MS scheme is straightforward 
in moment space: one simply divides uj(N,e) by <ft?-r^ (N,e), the square of the MS quark 




(2.23) 



F dis = ~^C F b (4 + C(2) - 2C(3)) + SFgl 



(2.24) 
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distribution, defined by 




^(«(e 2 ))+^(«(e 2 )) 



+ 0(1/N) . 



(3.1) 

with A the same function appearing in Eq. (2.15), while B$ is the virtual part of the non- 
singlet quark-quark splitting function (B$(a) = B(a) at one loop). As appropriate for 
an MS parton density, one can easily verify that 4>y^ (N, e) in Eq. (3.1) is a series of pure 
counterterms. Q is the factorization scale, which for simplicity throughout this paper we 
set equal to the Drell-Yan invariant mass 2 . We will now factor this density into virtual 
and real parts 

<j> m {N,€) = M^<l>R(N,€), (3.2) 



in such a way that in the (finite) ratio 



)_ (N) = 



( \r(Q 2 ,e)\ 2 

V Me) 2 



j; R (N,e) 2 U R (N,e) 
MN,e) 2 



+ 0{l/N) (3.3) 



the ratios of virtual functions and of real functions, displayed in the large brackets, are 
separately finite. To be precise, the factorization in Eq. (3.2) is uniquely defined by the 
following criteria: first, the ratio of virtual functions must be finite; second, as we are 
factorizing a series of pure counterterms, we would like also <f>v( e ) to consist only of poles. 
The real part 4>n(N,e) is then defined by Eq. (3.2). Note that 4>v(e) defined in this way 
is process-dependent, in contrast to cft^ (N, e); note also that, while (py^ has only simple 
poles of collinear origin, the real and virtual contributions will have cancelling double poles. 
We will now analyze separately the real and virtual contributions to the cross section. 

3.1 Cancellation of virtual poles 

The timelike Sudakov form factor has imaginary parts, which are the source of the largest 
contributions to the ratio in Eq. (2.20), while the MS distribution is real. We can thus 
simplify our analysis by writing 



|r(Q 2 ,6)| 
Me) 2 



r(Q 2 ,e) 



r(-Q 2 ,e) 



r(-Q 2 ,e) 

<Me) 



(3.4) 



The benefit of this lies in the fact that the first factor on the right hand side is finite, and 
already explicitly resummed in Eq. (2.20). The second factor, on the other hand, is purely 
real. Inspired by the explicit expression for the form factor, Eq. (2.19), we try the ansatz 



4>v{e) = exp 



1 



de 



~ 1 wx 2 

K{a s ,e) + G{a{e)) + l - / (a (A 2 )) 



(3.5) 



which has the same structure as the Sudakov form factor in Eq. (2.19), with the difference 
that G has no order e terms. 



2 It is straightforward to repeat the analysis below keeping these scales different. 
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We will now show, to all orders in the strong coupling, that the perturbative coefficients 
of the function G can be chosen so as to render Eq. (3.4) finite, and we will provide for them 
an explicit construction. Since the first ratio on the right hand side of Eq. (3.4) is finite, it 
is sufficient to prove the cancellation of poles for the second ratio, T(—Q 2 , e)/0y(e). First 
of all, one observes that all divergences arising from the terms K and manifestly cancel 
between r(— Q 2 ,e) and 0y(e), since they come from limits of integration independent of 
the particular energy scale considered. Hence divergences could only be generated by the 
G terms. To show that those terms can be made finite as well, we begin by writing 

oo oo 

G( as ,e)=J2J2 G ^ em {^y > ( 3 - 6 ) 

n=l m=0 

and by noting that the integration over the energy scale can be rewritten as one over the 
running coupling making use of 

dji 2 _ 2 da _ 1 da 1 

/i 2 P(a,e) e a 1 ^ , ( a\m ' 

1+ 4e ^ ° m -! {■*) 
m=l 

To proceed, let us define the truncated perturbative expansions for the various functions 
involved by 

n=l 

M oo 

G m (a., e) = E E G - £m ( v)" ' (3 ' 8) 

n=l m=0 

M 

o ( M,( a ,) S E<5 , " , (v)"- 



n=l 

Finally, define 

<(«.,«) s ^Mi^illz^aM , (3. 9 , 

P(M-i){ai s ,e) 

where the index (p) denotes the fact that the perturbative expansion for the ratio R is 
truncated at order a p s . 

Writing down the integrands of both the numerator and the denominator of the second 
ratio in Eq. (3.4) and keeping in mind the overall factor of 1/e from Eq. (3.7), one can 
easily formulate our goal in terms of R. We need to show that G can be chosen so that 

RM(a a ,e)=0(e) , VM . (3.10) 

The proof proceeds by induction. First of all, one sees immediately that it is possible to get 
R^(a s ,e) = 0(e), simply by choosing G^ = G^ . Next we assume that Eq. (3.10) holds 
for some fixed order M in perturbation theory for chosen values of , up to n = M, and 
show that G( M+1 ) may be chosen so that Eq. (3.10) holds also to order M + 1. 
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To isolate the terms of order <Jf +1 in -R^+jj , we use the identity 

1 1 P(M)(a s ,e) -d{M-i){ot s ,e) 



a,, e 



Substituting Eq. (3.11) into the definition of -R^+Jj, and neglecting all terms 0{a^ +2 ) 
and higher, as well as terms which are explicitly 0(e), one finds 



(3.11) 



R 



(M+l) 
(M+l) 



(a s ,e) = ^ +1) (a s ,e) + (gJ M+1) " - C^ 1 )) . (3.12) 



The theorem is now proven if one can show that R^^t^ has no poles in e, since a constant 

at 0(a^ +1 ) can be removed by suitably defining & M+l \ To see that this is the case, note 
that Ry^' 1 \a s , e) must, by the induction hypothesis, be of the form 



M oo 



m=l 1=1 



P=P0 



Multiplying this times /3(M-i)( a s, e), we must get back the numerator of R^-^K up to 
0(a^ +2 ) corrections, namely 



n=l 

M oo 



V 7T 



M oo 



.m=l Z=l 

M 



P=P0 



n=l m=0 n=l 



(3.14) 



Now, the right hand side has no poles in e, and no term of order a^f +1 , thus on the left hand 
side all poles and all terms of order a^ +1 must cancel. If po < 0, on the other hand, one 
will generate a term of the form d^i(a s /ir)( M+1 ^ (1/e) on the left hand side, which cannot 



be cancelled: in fact, note that the first sum in R 



(M+l) 
(M) 



starts at 0(e) by the induction 



hypothesis, and Pm-i has only a simple pole in e. We conclude that po > 0, as desired. 

To find explicit expressions one has to implement the observation that the coefficient 
of {(a s /ir) M+1 , e } in ^2^+^ must vanish. One verifies that G^ M+1 ^ is given by the 



coefficient of {(a s /ir 



e } in the ratio 



^(Af+i)(«s, e) = 



G(m+i) («s,e) 
A?(M)(«s,e) 



(3.15) 



Explicitly, 



g(M+l) = ^(M+l) 6q c (M) bl g (Af-l) + jo g (M-l) ^ G (M-2) + Ml G (M-2)^ ^ ^ 

o 4 ^ 4 ^ 16 2 4 ^ 8 2 
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3.2 Real emission contributions 

The complete expression for the MS -scheme DY cross section in the present framework is 
given by 



y AlS 



(TV) 



r(Q 2 ,e) 



r(-Q 2 



r(-Q 2 



Vie 



U R (N,e) 



We) 



(3.17) 



where the factor in square brackets arises from real gluon emission. Each function appearing 
in the real emission contribution exponentiates: ipR and Ur according to Eqs. (2.5) and 
(2.6), respectively, while <f>R is defined as the ratio of Eq. (3.1) to Eq. (3.5). Renormalization 
group arguments can be applied to each function, in d = 4 — 2e dimensions, as described 
in the Appendix. It is interesting to notice that running the coupling in d dimensions 
generates poles at two loops in the ratio UriPr/4>r, although the input at one loop is finite. 
The ratio must however be finite to all orders, as a consequence of the factorization theorem, 
together with the finiteness of the virtual contributions demonstrated above. This poses 
constraints on the two-loop coefficients of the functions involved, as described in more 
detail in Section 4, tying together real and virtual contributions to the cross section. 

Collecting and organizing the exponential expressions of the real emission functions, 
one can cast Eq. (3.17) in the standard form 



-'Ms 



(N) 



r(Q 2 ,e) 



r(-Q2 



x exp 



dz 



■ m-Q\e) 

\4>v(Q 2 ,e 

z N-l _ 1 



1 



+ D (a s ((1 - zfQ 2 )) jj +0(l/N) . 
A one loop calculation, with the inclusion of running coupling effects, yields 




(3.18) 



log 



r(-Q 5 



<MQV) 



IT 



C(2) 



a, 
I — 

7T 



C F b 



1 - — C(2) - —((3)] + 8RW 
32 sv ' 24 sv ' 



F MS ( a °) = —°F 



TV 



-\b C F (l - |c(2) - JC(3) 



+ SE. 



(2) 



MS 



(3.19) 



where 5R^ and are genuine two-loop contributions, formally unrelated to running 

coupling effects. Note that the function D in Eq. (3.18) is the same as in Eq. (2.15): 
a non-trivial statement, due to the fact that such a function, summarizing wide angle 
soft radiation, can be taken to vanish in the threshold-resummed deep-inelastic structure 
function [15, 25, 26, 27]. 

We can now turn to a discussion of the relevance of the exponentiation of constants in 
the two schemes, as tested at the two-loop level. 

4. Usage and limits of exponentiation 



It is clear that the results we have derived on the exponentiation of AMndependent terms 
do not have the predictive strength of the standard resummation of threshold logarithms. 
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In that case, in fact, the pattern of exponentiation is highly nontrivial, and the perturbative 
coefficients of entire classes of logarithms can be exactly predicted to all perturbative orders 
performing just a low order calculation. In the present case, even though constant terms 
exponentiate, the determination of the exact value of the AMndependent contribution at, 
say, g loops, always requires a g-loop calculation, albeit in some cases a simplified one. 

It remains true, however, that exponentiation and running coupling effects generate 
contributions to all orders which originate from low-order calculations. These contributions 
are an easily computable subset of higher order corrections, and it is reasonable to use them 
to estimate the full impact of higher orders. Specifically, given a (7-loop calculation of one 
of the cross sections we have discussed, one can extract the value of the various functions 
appearing in the exponent to that order, and then use exponentiation and RG running to 
estimate the (g + l)-4oop result. Given the existing results at two loops [28, 29], one could 
construct an estimate of the three-loop partonic cross section. Before embarking in such a 
calculation, it is however advisable to test the case g = 1, i.e. to use the two-loop results 
to verify the reliability of the method, by comparing the exact results with the estimate 
obtained by exponentiating the one-loop calculation and letting the couplings run. To this 
end, we can expand the partonic cross section in scheme s, u) s (N) as 



00 

(iV) = J>W(iV)(^) P . (4.1) 



p=0 

Next, we identify the coefficients of different powers of log AT, by writing 

2p 

4 P HN) = Y,^H^N + lE y . (4.2) 

i=0 

So far we are dealing with the exact cross sections. Let now uj'fj be the estimate for u^fj 
obtained by evaluating the exponent exactly at p — 1 loops, adding running coupling effects, 
and expanding the result to order p. One can define the deviation 

Jp) _ s (p) 

Acg - • (4.3) 

to ■ 

s,i 

In computing our estimates for the DIS scheme, we will employ Eq. (2.15), with the one- 
loop results taken from Eq. (2.17), Eq. (2.21) and Eq. (2.23). Furthermore, since we are 
taking into account running coupling effects, we will include the terms proportional to bo in 
Eqs. (2.24) and (A. 8). As far as the MS scheme is concerned, we will make use of Eq. (3.18), 
together with Eq. (3.19), with all purely two-loop contributions defined to vanish. Exact 
results are taken from Refs. [28, 29], focusing on "soft and virtual" contributions (all other 
contributions are suppressed by powers of N at large N) . To get numerical results we focus 
on SU (3) and set nj = 5. The results for the deviations Aw^' in the two schemes are given 
in Table 1. To gain a little further insight, we also separate the contributions proportional 
to the possible combinations of group invariants arising at two loops (i.e. C F , CaCf and 
nyCir), and display the results in Table 2 for the powers of log A" which do not lead to 
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i 





1 


2 


3 


4 


DIS 


0.26 


1.17 


0.13 








MS 


- 0.69 


1.79 


0.33 









Table 1: The deviations Aui s as denned in the text, for the DIS and MS schemes, 
exact agreement (i.e. i = 0,1,2). Here we give the coefficients separately for u^J and 

~(2) . ' 

lj s ■ for each scheme, since some of the exact coefficients vanish, so that the corresponding 
deviations as given in Eq. (4.3) are ill-defined. 





DIS 


MS 




estimate 


exact 


estimate 


exact 


coefficients of C\ 


i = 


38.06 


39.03 


3.33 


1.79 


i = 1 


- 13.09 


- 14.41 








i = 2 


9.85 


9.85 


5.16 


5.16 


coefficients of CaCf 


i = 


5.63 


18.12 


13.12 


6.82 


i = 1 


9.83 


- 0.25 


1.51 


- 0.47 


i = 2 


- 0.69 


0.35 





2.08 


coefficients of rijCp 


i = 


- 1.02 


- 4.40 


- 2.38 


- 0.80 


i = 1 


- 1.79 


0.35 


- 0.27 


- 0.52 


i = 2 


0.12 


- 0.15 





- 0.56 



Table 2: Comparison between estimates from exponentiation and exact results at two loops, 
presented separately for different colour structures, both in the DIS and MS schemes. 

Several remarks are in order. From Table 1, we observe that, as expected, leading 
(i = 4) and next-to-leading (i = 3) logs at two loops are exactly predicted by one-loop 
results and running coupling effects. Similarly, also as expected, NNL logarithms (in this 
context i = 2) have a small discrepancy which is entirely traceable to the two-loop cusp 
anomalous dimension 7^. At the level of iV-independent terms (i = 0), the agreement is 
quite reasonable, and in fact rather satisfactory in the DIS scheme, where exponentiation 
accounts for three quarters of the exact answer. 

Single log results, on the other hand, are much less satisfactory, displaying a discrep- 
ancy larger than 100% in both schemes. The reasons for this discrepancy are slightly 
different in the two schemes, but they highlight the same generic problem. In the DIS 
scheme, as discussed in the Appendix, running coupling effects in the ratio of parton dis- 
tributions ip/x generate a single log term with the wrong sign with respect to the exact 
result. This term must then be compensated by contributions which we would describe as 
'genuine two-loop', which as a result display a rough proportionality to bo. 

This phenomenon could be described as an 'excess of factorization', in the following 
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sense: to achieve the accuracy and generality of Eq. (2.15) it is necessary to introduce 
several functions, depending on different scales. Not all of these dependencies, however, 
are physical, and there may be (in fact there are) large cancellations in the scale dependence 
between different functions. This fact has been observed in the past: Catani and Trentadue 
pointed it out in Ref. [30], and more recently a similar observation was made by Gardi 
and Roberts in Ref. [27]. As a consequence, approximate coefficients that are dominated 
by running coupling effects, but not completely determined, may turn out to be quite 
inaccurate. 

In the MS scheme, the same kind of cancellation is displayed in a different way: there, 
as described in more detail below, one may use the constraint imposed by the finiteness 
of the real emission contributions to reexpress the single log coefficient in terms of purely 
virtual functions, and in the process the weight of running coupling effects changes con- 
siderably. Again, this indicates that computed running coupling effects may easily be 
compensated by unevaluated two-loop contributions. 

Finally, one may observe in Table 2 that abelian (C F ) contributions exponentiate with 
impressive accuracy, particularly in the DIS scheme. The slight superiority of the DIS 
scheme in this regard appears to be a fairly generic feature in Tables 1 and 2, perhaps to 
be ascribed to the more direct physical interpretation of the subtractions, as compared to 
the MS scheme. 

To conclude, we would like to point out that the methods of exponentiation we have 
outlined may be used not only for numerical estimates, but also to obtain, or test, analytical 
results. Specifically, since all functions employed have precise diagrammatic definitions, the 
computation of certain coefficients at two loops may be simplifed using this approach, as 
compared to a full calculation of the cross section. Further, the factorization into separately 
finite real and virtual contributions leads to constraints connecting different coefficients, so 
that different two-loop results can be nontrivially connected. To give an example, consider 
Eq. (3.17). There, the finiteness of the ratio Ur^r/4>r at two loops imposes constraints 
tying together real and virtual contributions to the cross section (recall that 4>R is defined 
as <t>yi$/4 > v)- Using the methods outlined in the Appendix and imposing the cancellation 
of double poles in the ratio of real functions, one may verify that the two-loop coefficient 
of log 2 N in loj^(N) must equal the two-loop cusp anomalous dimension 7^, as is well 
known. Further, imposing the cancellation of single poles in the same ratio, one finds 
that the value of the function D at two loops is completely determined by purely virtual 
diagrams. One finds 

£>( 2 ) = ^((2)b C F + 4E?f - , (4.4) 

(2) 

where B\ , the two-loop virtual part of the non-singlet quark-quark splitting function [31, 
32], is given by 

B f » = lei ( 1 - i« 2 ) + C(3,) + ^ (£ + >) - 3C<3>) - ( \ + «,)) , 

(4.5) 

while the second order contribution to the function G can be determined via Eqs. (3.16) 
and (2.19), by matching the resummed expression for T(Q 2 , e) to the explicit results for the 
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dimensionally regularized one- and two-loop Sudakov form factors, as given for example 
in Ref. [8, 33, 34]. We find 



& 2) = G (2) _ &0 G (1) 

<#> = 3 <?J± - km + C(3)l - ^ r 13 - " C (2) - ^#C(3) 



16 2 >w ^ v 7 4 V 3 ^ v y 108 

n f C F /209 . A 

'^T + C(2) , (4.6) 



6 V 36 



G? =C f (a- ^C(2)) 



We have thus rederived the coefficient D^ 2 \ obtained earlier in Refs.[13, 14] through match- 
ing to the two-loop cross sections of Ref. [28, 29], by using only information from purely 
virtual contributions. 



5. Conclusions 



We have shown how to organize all constants in the N— th moment of the Drell-Yan cross 
section in the DIS and MS schemes into exponential forms. Our MS -scheme result has 
the special feature that real and virtual contributions are separately finite. This organiza- 
tion rests crucially upon the refactorization properties of the unsubtracted Drell-Yan cross 
section and, for the DIS scheme, of the non-singlet deep-inelastic structure function, near 
threshold [3] . For the MS scheme the organization involves the construction of an exponen- 
tial series of pure counterterms that cancels all divergences in the spacelike Sudakov form 
factor. We have proven this cancellation to all orders. We emphasize that our arguments 
imply exponentiation to the same degree of accuracy for the MS -scheme DIS cross section 
]?2(N), although we have not given a detailed evaluation in that case. 

Although exponentiation of ^-independent terms does not have the same degree of 
predictive power as the resummation of threshold logarithms, it can be used with some 
degree of confidence to gauge the impact of higher order corrections to fixed order cross 
sections: we found that ^-independent contributions at two loops are reasonably well 
approximated by the exponentiation of one-loop results. The refactorization approach 
also leads to nontrivial connections between real and virtual contributions to the cross 
section, which can be used to test or in some cases simplify finite order calculations. On 
the negative side, one cannot in general trust running coupling effects to give by themselves 
a good approximation of two-loop results, unless the various scales at which the couplings 
are evaluated are tied to the physical scales of the full cross section. 

It might be interesting to make full use of the available two-loop information for the 
Drell-Yan and DIS cross sections to provide an estimate of three-loop effects, along the lines 
of Section 4. We prefer to regard the techniques presented here, which extend threshold 
resummations to a new class of terms, as a step towards the analysis of yet other classes of 
perturbative corrections which might be expected to exponentiate. A natural example is 
given by threshold logarithms suppressed by an extra power of the Mellin variable N, which 
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have recently been analized in the case of the longitudinal DIS structure function Fl(N) in 
Refs. [35, 36], and were shown to be phenomenologically important in Ref. [37]. Another 
possible extension of our work is of course the study of A-independent contributions to 
more complicated hard cross sections, involving more colored particles, along the lines of 
Refs. [38] and [39]: this would be an ingredient towards a precise resummed determination 
of such cross sections, which would be of considerable phenomenological interest for present 
and future hadron colliders. 
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Appendix 

Let us briefly discuss the application of renormalization group techniques to unconven- 
tional parton distributions such as the functions if)(N, e) and x{N, e) described in the text, 
specifically focusing on the contributions involving real gluon emission. The techniques of 
Ref. [3] lead to an expression of the form 

1>r(N, e) = exp y\ z ^—j'^- ^ ( ^-ti Q ,a a tf), e) } + O (1/N) , (A.l) 

and similarly for \- The functions i/jr and \R are both renormalization group invariant 
(i.e. their respective anomalous dimensions vanish) to this accuracy, for slightly different 
reasons: ipR cannot have overall UV divergences because its phase space is restricted to 
fixed total energy emitted in the final state. This automatically restrict also transverse 
momentum, so the phase space integration is UV finite. \Ri on the other hand, has a 
phase space restricted to fixed total light-cone momentum fraction, so that in principle 
it may have UV divergences arising from transverse momentum integrations. These di- 
vergences are in fact present, however it can be shown that, at least at one loop and in 
the chosen axial gauge, these divergences are suppressed by powers of A. Note that this 
is not in contradiction with the fact that the divergent terms for any quark distribution 
must be proportional to the Altarelli-Parisi kernel. It simply means that the corresponding 
divergences are of IR-collinear origin for the distributions at hand. 

The consequence of this statement for the functions and k x is that 

where (3 (e, ct s ) is the f3 function in d = 4 — 2e. An identical equation is obeyed by k x . Such 
equations can be solved perturbatively to determine the dependence of the distributions on 
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the momentum scale. Consider again for example n^, and let £ = (1 — y)Q/fj,. Expanding 



n=l 

one easily finds that at the two loop level the coefficients must be of the form 

,(1) fc A - J 1 ) (A <-2e 



4 ; «.o = <iwr 



4 2) & e ) = 4 2 i w r* + b £C w r 2e (r 2£ - 1) . f a . j ) 



explicit evaluation at one loop [3] yields 



K % (e) = 2C F (4nY IjL-± , 



k { £ (e) = 2C F (47r) e T(2 + e) cos (ire) . (A.5) 

As observed in Section 2, the finiteness of the ratio i^r/xr^ which is a consequence of fac- 
torization, requires that K^ Q (e) — ^^(e) = 0(e 2 ), since the double integration in Eq. (A.l) 
generates a double pole. In fact, upon redefing \i according to the MS prescription to absorb 
factors of log(47r) and je, 

*$> ( £ ) " *$ ( £ ) = 2C ^ [2 + C(2) + e (4 + C(2) - 2C(3)) + 0(e 2 )] . (A.6) 
One can go slightly further and observe that the finiteness of the ratio iPr/xr also constrains 

(2) 

the form of the pure two-loop contribution to and k x , given by the functions K^ Q (e) 

(2) 

and K x ' Q (e). Specifically, inserting Eqs. (A. 3) and (A. 4) into Eq. (A.l), and doing the same 
for xr, one finds that the ratio i/jr/xr will develop a simple pole in e at two loops, unless 



4,o (<o - «S ( £ ) = 2 Cpboe (2 + c(2)) + eH ^ + ° ( eS ) ' (A - 7) 

(2) 

in analogy with Eq. (A.6), with 5k 2 a constant arising at two loops to be used below. 
This constraint also fixes the coefficient of a contribution to the ratio proportional to log N 
at two loops, i.e. at NNL level. To be precise one finds 

3 „ , _„,A „ / 1 



exp 



^-C F (2 + C(2)) + (^) 2 (^4 2) + ^6 (2 + C(2)) (log AT + lE ) 



- ^C F b (A + ((2)-2C(3)))+OU,-,a 3 s 



(A.8) 



Once again, the contributions arising at two loops should be taken with a grain of salt when 
constructing the full cross section. It is true in fact that in this way we have determined the 
leading logarithmic contribution to this particular ratio, however in the full cross section 
there are competing log A?" terms arising at two loops from other functions, and in fact in 
the present case the logarithmic term in Eq. (A.8) goes in the wrong direction to 'predict' 
NNL logarithms at two loops, as discussed in the text. Similarly, there is no guarantee that 

(2) 

the uncalculated constant 5k 2 will not overwhelm the running coupling effects explicitly 
displayed in Eq. (A.8). 
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